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Abstract: There are two length-scales present simultaneously in all the principal directions for 
three-dimensionally expanding, finite systems. These are discussed in detail for the case of a 
longitudinally expanding system with a transverse flow and a transverse temperature profile. 
For systems with large geometrical sizes we find an mt- scaling for the parameters of the Bose- 
Einstein correlation function, which may be valid in the whole transverse mass region for certain 
model parameters. In this limit, the Bose-Einstein correlations view only a small part of the 
source. The large geometrical sizes can be inferred from a simultaneous analysis of the invariant 
momentum distribution and the Bose-Einstein correlation function. A preliminary analysis of 
the NA44 data indicates that instead of a small fireball we are observing a big and expanding 
snowball in S + Pb reactions at CERN SPS energies. 
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Bose-Einstein correlations and invariant momentum distributions are presented for ex- 
panding finite systems with applications to recent NA44 data. 

1. Introduction 

Bose-Einstein correlation functions (BECF-s) provide a unique tool for the analysis of 
the freeze-out geometry on the fermi scale. They are also subject to lots of non-ideal 
effects which make the analysis and interpretation of these data rather difficult. 

There are two length-scales present simultaneously in all the principal directions of three 
dimensionally expanding, finite systems |T|. One of the length-scales is the geometrical size 
Rg, the other is generated by the freeze-out temperature distribution and the gradients of 
the flow. This second kind of radius is referred to either as the 'thermal radius' Rt UHl 
or the 'length of homogeneity' The thermal radius characterizes the size of the 

region in the fluid from which particles with similar momentum are emitted, while the 
geometrical radius is present due to the finite size of the expanding system. 

2. Model specification 

We model the emission function for high energy heavy ion reactions as 

S(x, K) d x = 
u(x) ~ 
T(x) = 

Thus we include a finite duration, H{t) oc exp(— (r — t ) 2 /(2At 2 )). The decrease of 
the temperature in the transverse direction is controlled by the parameter a, while the 
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strength of the transverse flow is controlled by the parameter b. A non-relativistic limit 
for a = is described in ref. [[Ij]. One dimensionally expaning finite systems correspond 
to the a = b = case [fj], see this paper for the notation too. The first version of the 
three dimensionally expanding cylindrically symmetric finite model corresponds to the 
a = and b = 1 case [Q. The integrals of the emission function are evaluated using the 
saddle-point method The saddle-point equations are solved in the LCMS 0, the 

longitudinally comoving system, for 7] s « 1 and r xs << r . These approximations are 
warranted if | y — yo \« 1 + A7] 2 m t /T and f3 t = p t /m t « (a 2 + b 2 )/b. The flow is non- 
relativistic at the saddle-point if f3 t « (a 2 + b 2 )/b 2 . The radius parameters are evaluated 
here up to 0(r XtS /To) + 0(r] s ), keeping only the leading order terms in the LCMS. However, 
when evaluating the invariant momentum distribution (IMD), sub-leading terms coming 
from the cosh(r] — y) pre-factor are also summed up, since this factor influences the IMD 
in the lower m t region where the data are very accurate. 

2.1. Bose-Einstein correlations 

The BECF is parameterized in the form of C(Qi,Q sic ie,Qout) = 1 + A exp(— R\Q\ — 
RsideQlide ~ R-outQout) where the intercept parameter A and the radius parameters may 
depend on the rapidity and the transverse mass of the pair. We obtain |^]|| 

R 2 lde = Rl R 2 out = Rl + f3 2 Ar 2 , R 2 L = r 2 A V l (4) 

Rl = W T + W G A^ = A? + A^' (5) 

i.e. the parameters of the BECF are dominated by the smaller of the geometrical and the 
thermal length-scales. The thermal length-scales Rt and Ai]t are found to be 

t 2 T T 

Rt = 2 ° 72 — and a Vt = — • (6) 

a + b m t m t 

These analytic expressions indicate that the BECF views only part of an expanding source. 
Even a complete measurement of the parameters of the BECF as a function of the mean 
momentum of the particles may be insufficient to determine uniquely the underlying space- 
time emission function due to this reason. If the finite source sizes are large compared to 
the thermal length-scales and if we also have a 2 + b 2 ~ 1, one obtains an m t -scaling for 
the parameters of the BECF, 

9 9 9 9 T (a 2 + b 2 ) 1 

RLe - Rim - Rl -To—, valid for p t « - 2 (7) 

m t b b 

Note that this relation is independent of the particle type and has been seen in NA44 
data 0. This m^-scaling may be valid to arbitrarily large transverse masses with j3 t w 1 
if b 2 << 1. The finiteness of the expanding system reveals itself in an out-long cross term 
too which is next to leading order in the LCMS, being smaller than R 2 out — R 2 ide , which 
measured difference is very small M, alternatively explained by ref. M. 

2.2. Invariant momentum distributions 

The IMD plays a complementary role to the measured Bose-Einstein correlation function 
|]-[§. Namely, the width of the rapidity distribution at a given m t as well as T* the effective 
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temperature at a mid-rapidity yo shall be dominated by the longer of the thermal and 
geometrical length-scales. Thus a simultaneous analysis of the Bose-Einstein correlation 
function and the IMD may reveal information both on the temperature and flow profiles 
and on the geometrical sizes. E.g. the following relations hold: 

Ay(m t ) 2 = Ar] 2 + Ar) 2 T (mt), and 1 = __J + Izi. (8) 

i* 1 + lGKjrit = m) 1 

The geometrical contribution to the effective temperature is given by T G = T R 2 G jR\-. and 
the fraction / is defined as / = b 2 /(a 2 + b 2 ), satisfying < / < 1. The saddle-point sits 
at T) s = {y - y)/(l + Arf / Ar]^) , r X)S = (3 t bR 2 J {tqAt^) and r y , s = 0. 

For the considered model, the invariant momentum distribution can be calculated as 



d n 



dy dm t (2ir) 




exp I ^ I m t (27tA V 2 t 2 ) 1/2 (2tt^) cosh(^) exp(+A^/2) x 
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This IMD has a rich structure: it features both a rapidity-independent and a rapidity- 
dependent low-pt enhancement as well as a high-pt enhancement or decrease. 

Within this model, the rapidity-independent low-pt enhancement is a consequence of the 
transverse mass dependence of the effective volume, from wich particles with a given mo- 
mentum are emitted. This is described by the (27tA?7^Tq ) 1 ^ 2 {2t[R 2 ) cosh(^ s ) exp(+A^/2) 
factor, being proportional to (T/m t )^ 3 ^ 2 ^ exp(+T/(2m t )). 

The rapidity- dependent low-pt enhancement, which is a generic property of the lon- 
gitudinally expanding finite systems 0, reveals itself in the rapidity-dependence of the 
effective temperature, defined as the slope of the exponential factors in the IMD in the 
low-pt limit at a given value of the rapidity. The leading order result is 

T TT ( T\ 2 

T *" M = i + „(»-»,)' with ° = w( A " 2 + ^ • (10) 

The high-pt enhancement or decrease refers to the change of the effective temperature 
at mid-rapidity with increasing m t . The large transverse mass limit shall be in general 
different from the effective temperature at low pt given by T* since 

rp 2T ^ 2 / T G (m) \ 

2- f % 2-f\ T + T G (m) J 

Utilizing Tq/T = R 2 G jR\, the high-pt enhancement or decrease turns out to be controlled 
by the ratio of the thermal radius Rrimt = m) to the geometrical radius Rq. One obtains 
> T* if Rr(m) > R% and similarly Too < T* if R^im) < R 2 G . Since for large colliding 
nuclei R G is expected to increase, the high-pt decrease in these reactions becomes a 
geometrical effect, a consequence of the large size. 
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3. Large Halo: Applications to NA44 data 

The effects of long-lived resonances on both the IMD and the BECF can be taken 
into account analytically along the lines of [IIJ, supposing that the halo is sufficiently 



large. NA44 data for central S + Pb reactions at CERN SPS with 200 AGeV show an 
approximately m t independent intercept parameter: X n + = 0.56 ± 0.02 and 0.55 ± 0.02 
at m t = 150 MeV and 450 MeV, respectively. This suggests that the halo contains 
1 - V\ = 25 ± 2 % of all the pions. If A = const the equations simplify JH| as 

1 - r ^ = -r-r2 = -r-T2> and C(Ak,K) = i + \R c (Ak,K). (12) 

ay am t ay am t ay am t 

In this case the only apparent effect of the halo is to reduce the intercept parameter of 
the measured BECF to A = const while the IMD and the (Ak, K) dependence of the 
BECF is determined by the central part apparently exclusively. This central part is well 
accessible to the analytical models, presented e.g. in the previous sections. In general the 
halo may have more influence both on the IMD and the BECF. 

In a preliminary analysis, we have fitted the NA44 preliminary IMD for pions and kaons 
together with the final NA44 data for the mt dependence of the BECF parameters for 
both pions and kaons. Fixing the parameters yo = 3, a = and 6=1 and extrapolating 
the equations to the relativistic m t domain we get a description of the IMD and the BECF 
for both pions and kaons at a x 2 /NDF = 2.2. The preliminary analysis indicates large 
geometrical source sizes for the pions, Rg{^) ~ 8 fm, late freeze-out times of To ~ 8 fm/c, 
a vanishing duration At for both pions and kaons and finally surprisingly low freeze-out 
temperatures, T « 80 MeV. The kaons appear from a smaller region of Rg{K) « 3.5 fm. 



4. Conclusions: 

Instead of observing a small fireball we are observing a big and expanding snowball at 
CERN SPS S + Pb reactions accoringly to the above preliminary analysis of the partly 
preliminary NA44 data. The snowball may contain a hot core, to be investigated in a 
future data analysis. The preliminary results are very sensitive to the detailed structure 
of the m r dependence of the parameters of the Bose-Einstein correlation functions as well 
as to the details of the invariant momentum distribution of both pions and kaons. 
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